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We derive bounds on the dark matter annihilation cross-section for low-mass (5 — 20 GeV) dark 
matter annihilating primarily to up or down quarks, using the Fermi-LAT bound on gamma-rays 
from Milky Way satellites. For models in which dark matter-Standard Model interactions are 
mediated by particular contact operators, we show that these bounds can be directly translated into 
bounds on the dark matter-proton scattering cross-section. For isospin-violating dark matter, these 
constraints are tight enough to begin to constrain the parameter-space consistent with experimental 
signals of low-mass dark matter. We discuss possible models that can evade these bounds. 
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Introduction. There has been great interest in low- 
mass dark matter (mj ^5 — 20 GeV) as a possible 
way to explain potential signals from the DAMA pQ , Co- 
GeNT [2] and CRESST [3] experiments. However, this 
data seems to be in tension with exclusion bounds from 
XENON10 |4], XENON100 [5], CDMS [7] and SIM- 
PLE [8]. Of these, the constraints from the xenon-based 
experiments appear to be the tightest. Attention has 
thus turned to new models which can potentially recon- 
cile this data. One focus has been on isospin-violating 
dark matter (IVDM) [9-15J, in which dark matter in- 
teractions between protons and neutrons are different. 
Although further reconciliation is required for complete 
consistency |16) . destructive interference between dark 
matter couplings to protons and to neutrons can poten- 
tially relieve the tension between the CoGeNT, DAMA 
and XENON10/100 data sets. IVDM may thus play an 
important role in understanding the low-mass data. 

A consequence of destructive interference is that the 
coupling of dark matter to up- and down-quarks must 
be relatively large for a given value of the dark matter- 
nucleus spin-independent scattering cross-section. This 
translates into larger annihilation and production cross- 
sections. Indirect detection and collider searches, where 
destructive interference plays no role, can thus provide 
stronger constraints on such dark matter models. In par- 
ticular, the large annihilation cross-sections can enhance 
photon signals from the decays of hadrons produced in 
dark matter annihilation to up- and down-quarks. Be- 
cause of their sensitivity to dark matter with mass in 
the range ^5 — 20 GeV, gamma-ray searches of dwarf 
spheroidal galaxies [TThTT?] can constrain these models. 

In this Letter, we determine the constraints on the dark 
matter annihilation cross-section arising from gamma-ray 
searches of dwarf spheroidals, assuming dark matter an- 
nihilates entirely to up or down quarks. We then apply 
these constraints to IVDM models. We will see that if 
isospin- violation is the key feature in alleviating the ten- 
sion between the DAMA, CoGeNT and XENON10/100 



data sets, then models with purely contact interactions 
are tightly constrained by Fermi-LAT data. We also dis- 
cuss how constraints on IVDM models weaken when dark 
matter-quark interactions occur through a light media- 
tor. 

Bounds from dwarf spheroidals. In [18] , bounds 
on the dark matter annihilation cross-section were de- 
termined from a combined analysis of gamma-rays from 
dwarf spheroidals. This analysis compared the observed 
number of high-energy photons arriving directly from the 
dwarf spheroidals to the observed background slightly 
off-axis from the dwarf spheroidal. These bounds were 
expressed in terms of the quantity $pp, defined as 




where Bf is the branching ratio for dark matter to an- 
nihilate to channel /. dNf/dE is the photon spectrum 
for each annihilation channel and E t hr is the detector's 
threshold photon energy. Using Ethr = 1 GeV, [TS] 
found the 95% CL bound 

$ PP < 5.0±4'j x 1(T 30 cm 3 s" 1 GeV" 2 , (2) 

where the asymmetric uncertainties are 95% CL system- 
atic errors |17j . resulting from uncertainty in the mass 
density profile of the various satellites. 95% CL exclusion 
contours on the dark matter annihilation cross-section 
can thus be determined from eqs. |1|2| for any mj, using 
the integrated photon distribution from each annihilation 
channel. 

For annihilation to up or down quarks, we generated 
the photon spectrum using Pythia 6.409 [20], run on the 
computing cluster at the Hawaii Open Supercomputing 
Center. The up and down quark channels produce more 
photons than the r-channel, resulting in tighter bounds. 
The bound for the uu channel is plotted in figure [lj the 
bound for the dd channel is slightly tighter than that for 
the uu channel, but the difference is never more than a 
percent for the given mass range. For these bounds, we 
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FIG. 1. 95% CL exclusion bound (see text) from Fermi in the 
(mx, {uav}) plane, assuming dark matter annihilation to uu. 
The results for annihilation to the dd channel are visually 
identical. Also shown is the exclusion contour from BESS- 
Polar II for the same channel (assuming a diffusion halo size 
of 4 kpc). The dotted line indicates 1 pb, which is approx- 
imately the total annihilation cross-section at freeze-out for 
dark matter produced thermally in the early universe. 



have assumed the central value of $pp < 5.0, applica- 
ble for the satellite mass density profiles found in [17] : 
allowing this to vary within 95% systematic error bars 
could weaken this bound by a factor of up to ~ 1.9, or 
strengthen it by a factor of up to ~ 10. 

Dark matter annihilation to up- and down-quarks can 
also be constrained by limits on the cosmic ray anti- 
proton flux from experiments such as BESS-Polar II [21] 
and Pamela [2"2"] 123 ] . Of these two experiments, BESS- 
Polar II reports tighter bounds, which we have plotted 
in figure [I] (light blue dashed curve, assuming a diffu- 
sion halo size of 4 kpc). The reported BESS-Polar II 
bound on dark matter annihilation to up- or down-quarks 
is tighter than the Fermi-LAT bound, given the astro- 
physics assumptions underlying both analyses. However, 
there are significant systematic astrophysical uncertain- 
ties which can weaken the anti-proton flux limits by up 
to a factor ~ 50 [24] • These include uncertainties regard- 
ing the background and the dark matter density profile 
of the galaxy, as well as uncertainties in the cosmic ray 
propagation model. Solar modulation can also sensitively 
affect bounds on dark matter annihilation arising from 
anti-proton flux measurements. Depending on these sys- 
tematic uncertainties, the constraints on dark matter an- 
nihilation from the anti-proton flux may be weaker than 
those from gamma-ray searches of dwarf spheroidals. 

The anti-proton flux uncertainties which we have dis- 
cussed above are largely independent of the systematic 
uncertainties for gamma-rays from dwarf spheroidals, 
which derive primarily from the dark matter distribu- 
tion within the spheroidal. These two indirect detection 
methods are thus quite complementary. 

IVDM. We now apply these constraints to models of 
isospin- violating dark matter. We will consider the case 
in which dark matter interactions with Standard Model 
particles (both scattering and annihilation) are mediated 



by effective contact operators. The most general such 
effective operator can be written in terms of dark mat- 
ter and quark bilinears, each of which transforms as ei- 
ther a scalar, pseudoscalar, vector, pseudovector, ten- 
sor or pseudotensor. For such interactions, the dark 
matter-proton spin-independent scattering cross-section 
(cgj) can be unambiguously related to the dark matter 
annihilation cross-section. 

We will be interested in IVDM models which couple 
only to up and down quarks, since these provide conser- 
vative bounds. Couplings to all other quarks are isospin- 
invariant, and would result in a larger dark matter an- 
nihilation cross-section for fixed <7gj. We will only con- 
sider operators which yield spin-independent, velocity- 
independent scattering matrix elements. Furthermore, 
one will only obtain relevant bounds from Fermi data if 
the annihilation matrix element is s-wave. Given these 
choices, there are only three effective operators to con- 
sider: 
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where Od,c,r axe possible couplings if the dark matter 
(X) is a Dirac fermion, complex scalar or real scalar, 
respectively^ Here, C q u c R are dimensionless couplings, 
M* is an overall energy scale and q — u,d. The WIMPless 
model [2"f5H2"B] exhibited in [13] as an example of IVDM 
was a real scalar coupling through Or. 
The fjgj™ can then be written as 
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where fi p is the dark matter-proton reduced mass and 
the f Pin are coefficients parameterizing the coupling of 
dark matter to protons and neutrons, respectively. For 

fn/fp 0.7, the constraints from XENON10/100 are 

maximally suppressed because of almost exact cance- 
lation between proton and neutron interactions in the 
xenon target [THUS]. Moreover, for this choice of f n /f P , 



1 Majorana fermions are not considered because there is no cor- 
responding operator producing both s-wave annihilation and 
velocity-independent scattering. Also, it is useful to note that the 
most popular Majorana candidate, the neutralino in supersym- 
metry, is at least marginally inconsistent with the observations 
from CoGeNT, DAMA, and SIMPLE [25]. 
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the regions favored by DAMA and CoGeNT are brought 
into alignment. 

The f pn are related to the coefficients C u,d by 
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f P = C m BP 
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(5) 



where B^f\ are the integrated nuclear form- factors [25] . 
When the dark matter is a scalar it couples through op- 
erator Oq,r to the scalar quark current, and the form 



factors are roughly given by B%, ' — B^" 01 "' ~ 6 and 
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4. When the dark matter is a 



Dirac fermion it couples through operator Od to the vec- 
tor quark current, and the appropriate form factors are 

r,p(vec.) 



, = 2 and B u 



gp(vec.) _ ^ 



The annihilation cross-section arising from each effec- 
tive operator can be written as 
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If the dark matter is a Dirac fermion or complex scalar, 
we will assume that the dark matter particle and anti- 
particle number densities in Milky Way satellites are 
equal. 

For any choice of f n / f p and choice of the interaction 
operator, one can use eqns. |4|5|6 to obtain constraints on 



Si 



from the bounds in figure |l| These constraints are 
shown in figure[2]for f n /f P ~ —0.7, assuming dark matter 
interacting through Oq (if dark matter couples through 
Or ot Od, the bounds would be tighter by a factor of 
2 or ~ 4, respectively). Also shown are the signal re- 
gions of DAMA [1] (3a, assuming no channeling [301 l3"T] L 
CoGeNT (90% CL) and CRESST 0] (2a), along 
with 90% CL exclusion contours from XENON10 g], 
XENON100 0, CDMS 0(7] and SIMPLE 0. As we 
can see from this figure, the CoGeNT region of interest 
would be excluded by the Fermi data, under the assump- 
tions made here. 

The CoGeNT experiment has recently reported a pre- 
liminary analysis indicating that their experiment may 
have more surface area contamination than originally 
thought [32l E3] ■ The effect of this correction would be 
to shift the CoGeNT region of interest to slightly larger 
mass and slightly smaller <7gj, though the magnitude of 
the shift can only be determined by a complete analysis 
from the CoGeNT collaboration. This shift would serve 
to weaken constraints from Milky Way satellites. 

Uncertainties from particle physics and astro- 
physics. There are a several ways in which the bounds 
from Fermi data can be weakened. For example, if the ef- 
fective interaction operator is generated by the exchange 
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FIG. 2. Favored regions and exclusion contours in the 
(ra x , o p ) plane for IVDM with f„/f p = -0.7. The solid black 
line is the bound from this analysis at 95% CL, assuming cen- 
tral values for the satellite density profile. 



of a mediating particle with mass ~ 1 GeV, then 
the momentum transfer during scattering interactions at 
any direct detection experiment (O(10 keV)) would be 
much smaller than the mediator mass [34 . As a re- 
sult, scattering interactions would still be described by 
an effective contact operator. However, the annihilation 
cross-section would receive an additional suppression of 
(m 4> /2m x ) 4: ~ 1Q- 4 - 10~ 5 . In this case, models with 
CTgj which could match the CoGeNT data would be un- 
constrained by limits from dwarf spheroidals by ~ 3 — 4 
orders of magnitude. 

In addition, bounds from Fermi would be insignificant 
if dark matter interactions were mediated by an opera- 
tor yielding velocity-independent spin-independent scat- 
tering, but p-wave suppressed annihilation. An exam- 
ple of such an operator would be (l/M 2 )XXqq, in the 
case where the dark matter is a fermion. Another ex- 
ample would be (l/M 2 )X*d (1 Xq^ fl q, if the dark mat- 
ter is a complex scalar. Dark matter coupling through 
these operators would have an annihilation cross-section 
which is suppressed by v 2 ; the constraints from Fermi on 
Ogj would thus be suppressed by more than 6 orders of 
magnitude. These models could potentially explain the 
low-mass data of DAMA, CoGeNT and XENON10/100, 
but would be unconstrained by the Fermi-LAT dwarf 
spheroidal search. 

Both p-wave suppressed annihilation and light media- 
tors would similarly weaken the bounds from anti-proton 
flux. In either case, the constraints from the anti-proton 
flux are suppressed by the same factor as the gamma- 
ray bounds, as discussed above. Such models matching 
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the CoGeNT, DAMA and XENON10/100 data would be 
unconstrained by bounds from BESS-Polar II. 

The dark matter-proton scattering cross-section 
needed to match low-mass direct detection data can be 
shifted by other particle physics or astrophysics uncer- 
tainties, such as uncertainties in the nuclear form factor 
or in the dark matter velocity distribution near the earth. 
To normalize with all of the observations, here we have 
used the standard halo model to calculate the WIMP ve- 
locity distribution (e.g. 35J); variations from this model 
may affect the results presented, especially for high mass 
targets [36] . 

Complementary searches. There are interesting 
constraints on the annihilation of low-mass dark mat- 
ter from WMAP-7 data [371 [38]. These constraints arise 
from the effect on standard recombination of the pho- 
tons produced by dark matter annihilation. In the mass 
range we consider (assuming w/d-channel annihilation), 
the bounds which we obtain from dwarf spheroidals are 
a factor of ~ 5 tighter than those arising from WMAP-7 
data. Note there are some systematic uncertainties in the 
WMAP-7 bounds regarding the effect of proton produc- 
tion on heating of the cosmic medium [35] ■ Again, these 
uncertainties are different from and independent of the 
uncertainties in dwarf spheroidal searches. Planck data 
is expected to significantly improve these CMB bounds. 

There have been several recent studies of constraints 
on low-mass dark matter arising from the neutrino flux 
yielded by annihilation in the sun [11] I39H43] . But neu- 
trino detectors cannot constrain dark matter which an- 
nihilates primarily to first generation quarks. The final 
state quarks hadronize, and these light hadrons typically 
lose energy and stop in the sun before decaying. They 
thus produce a neutrino spectrum far too soft to be dis- 
tinguished from the atmospheric neutrino background. 

The Fermi-LAT gamma-ray bounds are also com- 
plementary to recent constraints from collider monojet 
searches [44TI48] on effective contact operators coupling 
dark matter to quarks. However, applying previous re- 
sults to the case considered here is difficult, since they 
have typically assumed a dark matter coupling to all 
quarks. In particular, bounds in previous studies depend 
strongly on the existence of heavy quark operators, and 
are weakened considerably when only couplings to first 
generation quarks are present. 

To address this issue, we provide simple bounds for Ogj 
given couplings only to up and down quarks with /„/ f p — 
—0.7. These bounds are based on an ATLAS monojet 
search using 1 fb _1 of data [33] ■ We present bounds for 
two possible selection cuts on j£ T and on the pr of the 
leading jet: p T > 120 GeV, $ T > 120 GeV and p T > 
350 GeV, $ T > 300 GeV. In addition to these selection 
cuts, the ATLAS study also used an intermediate cut of 
$ T > 250 GeV, p T > 220 GeV. The bounds provided 
by this intermediate cut lie between those of the other 
two, and are not shown. For each choice of selection 
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TABLE I. Upper bounds on <Tgj in pb from ATLAS monojet 
searches, assuming the listed cuts on the leading jet pt and 
on the missing transverse energy. Masses are listed in GeV. 
The columns correspond to Dirac fermions, complex scalars, 
and real scalars respectively. 



cuts, the ATLAS collaboration determined a bound on 
the pp — > XX j cross-section attributable to new physics, 
based on a comparison of the number of observed events 
to the expected number of background events. 

For IVDM models, we generated the production cross- 
section for each of the effective operators in eqn. [3] us- 
ing MadGraph/MadEvent 5.1.3 50 . The selection cuts 
were imposed at the parton level. According to the anal- 
ysis of (47] |48j [51], approximately 40% of the events 
satisfying the cuts at parton level will satisfy the cuts 
once hadronization and detector effects are included. 
Accounting for this 40% selection efficiency, wc deter- 
mined upper bounds on C^'^/Af* for the scalar case and 

C U jflMl for the Dirac case. These translate directly 
into the upper bounds on <7gj given in Table [i] for the 
two sets of missing energy and pt cuts. Equivalent Teva- 
tron bounds from CDF monojet searches [52] are weaker 
than the ATLAS bounds (assuming the pt > 120 GeV, 
$ T > 120 GeV selection cut) by a factor of ~ 4. If dark 
matter is a real or complex scalar the monojet search 
sensitivity is significantly weaker than the sensitivity of 
Fermi-LAT and direct detection experiments. Only in 
the case of Dirac dark matter are collider bounds com- 
parable to the Fermi-LAT constraint. 

However, in the case of light mediators, collider bounds 
become significantly weaker. Assuming that the qq — > 
XX process proceeds through the s-channel, the collider 



cross-section scales as (p^ — m|) . For light mediators 
we have <C \p<j,\ > 2mx- This leads to a collider 
cross-section which is suppressed by at least (2mx) 4 , 
instead of vnT^ . The corresponding rescaling of the direct 
detection bound relative to the contact operator case is at 
least (2mx/m < j > ) 4: , weakening the bound significantly [48] . 
If ~ 1 GeV, the bound on <7gj would be raised by a 
factor of ~ 10 4 — 10 5 , placing relevant direct detection 
cross-sections for any species of dark matter well out of 
range of collider studies even at the LHC, at least in the 
conservative case of couplings only to light quarks and 
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through low-mass mediators 

Conclusions. We have shown that Fermi-LAT 
gamma-ray searches of Milky Way satellites can constrain 
the cross-section for low-mass dark matter to annihilate 
to up or down-quarks. These constraints are tighter than 
for models with annihilation primarily to b-quarks or 
r-leptons, due to the larger number of photons arising 
from annihilation to first-generation quarks. This is es- 
pecially interesting in the case of isospin-violating dark 
matter, where destructive interference between up and 
down quark couplings imply an enhancement in the an- 
nihilation cross-section for models which can match the 
low-mass data of DAMA, CoGeNT and XENON10/100. 
In particular, IVDM models whose effective contact inter- 
actions yield s-wave annihilation arc tightly constrained. 
Nevertheless, there arc IVDM models which can evade 
these Fermi-LAT bounds and other constraints from indi- 
rect detection and collider searches. These IVDM models 
involve either low-mass mediators, or effective interaction 
operators which yield p-wave suppressed annihilation. 

It is worth noting that microscopic dark matter models 
may involve interactions mediated by more than one ef- 
fective operator. To understand constraints on such mod- 
els from gamma-ray searches, one would need to know 
how the coefficients of the various interaction effective 
operators are correlated in any particular model. 

In the case where dark matter annihilates to first gen- 
eration quarks, the sensitivity of indirect search strate- 
gies depends on the spectrum of first generation quarks. 
The photon spectrum arising from the hadronization and 
decay of first generation quarks is peaked at low val- 
ues of Ery I mx ■ Future searches using Milky Way satel- 
lites would thus greatly benefit in any reduction in the 
gamma-ray analysis threshold, E t h r . 
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